arXiv:1504.04126vl [hep-ph] 16 Apr 2015 


IPMU 15-0049 


Observable dark radiation 
from cosmologically safe QCD axion 


Masahiro Kawasaki 0,6 , Masaki Yamada 0,6 , and Tsutomu T. Yanagida 6 

a Institute for Cosmic Ray Research, The University of Tokyo, Kashiwa, Chiba 277-8582, 

Japan 

b Kavli IPMU (WPI), TODIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 


Abstract 

We propose a QCD axion model that avoids the cosmological domain wall problem, 
introducing a global SU(3) f family symmetry to which we embed the unwanted PQ 
discrete symmetry. The spontaneous breaking of SU(3) f and PQ symmetry predicts 
eight NG bosons as well as axion, all of which contribute to dark radiation in the 
Universe. The derivation from the standard model prediction of dark radiation can 
be observed by future observations of CMB fluctuations. Our model also predicts a 
sizable exotic kaon decay rate, which is marginally consistent with the present collider 
data and would be tested by future collider experiments. 


1 Introduction 


After the discovery of the 126 GeV Higgs boson at the LHC [1, 2], we obtain the complete 
parameter set of the Standard Model (SM) of particle physics except for the remaining 
parameter, the strong CP phase [3, 4]. We only know the upper bound on the strong CP 
phase, |#qcd| ~ 10 -(10-11 ), which comes from the experimental upper bound on the neutron 
electric dipole moment [5]. This is unnaturally small, and thus we face a severe fine tuning 
problem called the strong CP problem. In addition to that, the observations of neutrino 
oscillations require new physics beyond the SM, and the origins of baryon asymmetry and 
dark matter (DM) are also long-standing challenges in cosmology and particle physics. One 
of the simplest solutions to those problems is a QCD axion model [6, 7, 8, 9] with right- 
handed neutrinos [10, 11]. The strong CP phase is cancelled by the vacuum expectation 
value (VEV) of axion [12], which is the pseudo-NG boson associated with the spontaneous 
symmetry breaking (SSB) of Peccei-Quinn (PQ) symmetry. The axion can be generated as 
a coherent oscillation in the early Universe and is a good candidate of cold DM [13, 14, 15]. 
Introducing heavy right-handed neutrinos, we can explain small neutrino masses via the 
seesaw mechanism. The out of equilibrium decay of the right-handed neutrinos can generate 
lepton asymmetry, which is then converted to baryon asymmetry via the sphaleron effect [16]. 

The above simple scenario may confront the problem of sizable isocurvature fluctua¬ 
tions [17, 18, 19] or the cosmological domain wall problem [20, 21]. If the PQ symmetry 
is spontaneously broken before inflation, it predicts a sizable isocurvature density perturba¬ 
tions due to quantum fluctuations of the axion field during inflation. The resulting amount 
of isocurvature perturbations is inconsistent with the observation of CMB fluctuations unless 
the energy scale of inflation is relatively small to suppress the amplitude of axion fluctua¬ 
tions. On the other hand, it is possible that the PQ symmetry is broken after inflation. In 
this case, the SSB of PQ symmetry results in formation of cosmic strings, and then at the 
QCD phase transition domain walls whose boundaries are cosmic strings are formed [22], 1 
The domain walls are stable and soon dominate the Universe, so that they spoil the success 
of the standard cosmological scenario. Although this is not the case if the domain wall num¬ 
ber is equal to one, it restricts axion models. In the DFSZ axion model [9], in particular, 
the domain wall number is equal to 6, so that it confronts the domain wall problem. 

In Ref. [24], Lazarides and Shah have proposed a novel mechanism that solves the domain 
wall problem. They introduced an additional continuous symmetry to which the unwanted 

^ven if the PQ symmetry is broken during inflation, the PQ symmetry may be restored due to the 
thermal effects after the reheating for sufficiently high reheating temperature suggested from the successful 
leptogenesis [23]. 
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discrete symmetry is embedded as a center of the continuous group. Since the degenerate 
vacua is connected with each other through the continuous symmetry, there is no domain 
wall problem in this case. 

In this paper, we consider a variant of the DFSZ axion model and introduce a global 
S77(3)f flavour symmetry for three families of quarks and leptons to realize the Lazarides- 
Shafi mechanism. The unwanted discrete PQ symmetry is embedded to the SU (3)f symmetry 
so that the vacua are continuously connected to avoid the formation of stable domain walls. 
Both S77(3)f and PQ symmetry are spontaneously broken at the same time, so that the low 
energy effective theory contains eight NG bosons, called fa /mi Ions [25, 26, 27, 28], as well as the 
axion. We find that they are thermalized and then decoupled from the thermal plasma after 
the SSB, and give a sizable contribution to dark radiation in the subsequent cosmological 
history [29, 30]. We find that the resulting amount of the dark radiation is consistent with the 
present observed value of effective neutrino number, and the discrepancy from the standard 
model prediction is measurable by future observations of CMB fluctuations. In addition, 
the SSB of the SU( 3)f symmetry predicts collider signatures in terms of flavour changing 
processes. In particular, the exotic kaon decay rate is marginally consistent with the present 
constraint and would be measured by future collider experiments. 

In the next section, we introduce a variant of the DFSZ axion model with a global 
SU( 3)f flavour symmetry for three quarks and leptons, and briefly explain the Lazarides- 
Shafi mechanism. In Sec. 3, we consider constraints and predictions of our model: cold axion 
DM, dark radiation, exotic kaon decay, and baryon asymmetry. Section 4 is devoted to the 
conclusion. 


2 Model 


We consider a variant of the DFSZ model with global SU (3)f x U (1)pq symmetry. The charge 
assignment is shown in Table 1, where y is the held responsible for the SSB of SU( 3)f x 
I/(1)pq. The right-handed neutrinos u c are introduced to realize the seesaw mechanism [10] 
and thermal leptogenesis [11], which we discuss in Sec. 3.4. Note that the SU( 3)f symmetry 
has anomaly, so that it has to be a global symmetry. The SU( 3)f x D(1)pq symmetry is 
spontaneously broken completely by the VEV of y, which we denote as 


(x) lj 


1 

7 ! 


diag (vi,v 2 ,v 3 ). 


( 1 ) 


The SSB predicts eight NG bosons called familons [25, 26, 27, 28] as well as pseudo-NG 
boson called axion [12] in the low energy effective theory. 
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Q 

u c 

d c 

L 

e c 

v c 

H u 

H d 

X 

SU (3) f 

3 

3 

3 

3 

3 

3 

6* 

6* 

6* 

U( 1)PQ 

1/4 

1/4 

1/4 

1/4 

1/4 

1/4 

-1/2 

-1/2 

1 


Table 1: Charge assignment for matter fields. 


The axion couples with the gluon held through the anomalous interaction: 
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where g s is the QCD gauge coupling constant, G afiu is the gluon held strength, and G c f ] u is 
its dual. The axion decay constant F a is given as 

1 

_ 

F 2 _ 

a 

In this paper, we consider the case that the 577(3)f x I/(1)pq symmetry is spontaneously 
broken after inflation. Since the hrst homotopy group is given by Hi [577(3)f x U(1)pq] = Z, 
cosmic strings form at the time of the SSB of the SU (3)f x U (1)pq symmetry [22], Then, after 
the QCD phase transition, the axion held effectively acquires a periodic potential through 
the QCD instanton effect. The non-perturbative QCD effects associated with instantons 
break I/(1)pq to in our model [21, 31, 32], so that domain walls may appear as bounded 
by the cosmic strings at the QCD phase transition. 

We define the domain wall number N dw as the number of disconnected vacua. Each 
cosmic string becomes the boundary of IVdw domain walls after the QCD phase transition. 
If the SU( 3)f symmetry was absent, it would be given by IVdw = 3 because the non- 
perturbative QCD effects break U{ 1)pq to Z A domain wall system with ATdw > 2 is 
stable and thus eventually dominate the Universe, in which case the subsequent cosmological 
history is highly inconsistent with the observations [20, 21]. However, we introduce the 
SU( 3)f symmetry to embed the Z^ symmetry into the continuous symmetry, so that the 
degenerate vacua are continuously connected with each other. 2 This implies that there is only 
one minimum along the axion held, so that each cosmic string is attached by only one domain 
wall (IVdw = 1) after the axion acquires the effective mass at the QCD phase transition. In 
this case, the domain wall is short lived due to their tension, and the subsequent cosmological 
scenario can be consistent with the standard cosmology [22, 24], 

2 In Ref. [33] they consider a GUT theory with S'O(IO) x SU(3)t x Z7(1)pq to realize the Lazarides-Shafi 
mechanism. They use the center of S'O(IO) x SU (3)f to embed the discrete symmetry into the continuous 
group. In this model, however, the PQ symmetry breaking as well as the SU (3)f symmetry breaking occur 
at the GUT scale. 
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Here, let us check that the symmetry is actually embedded to the SU(3){ symme¬ 
try. The {7(1)pq symmetry commutes with SU(3){, so that its discrete subgroup Z also 
commutes with it. Therefore, what we should check is that the Z^ group is embedded to 
the center of SU( 3)f because the center of a group commutes with all the group elements by 
definition. The center of S77(3)f group is Z 3 and its effect on matter fields can be written as 

Q ->■ e 2ni{n/3) Q 
H u ->■ e 2ni{ ~ 2n/3) H u 

^ _s. e i27r(—2rx/3) ^ _ e i2n(4n/3 )^ 

where n (= 1,2,3) is an integer. The other matter fields are also transformed in a similar 
way. Here, you can see that the center of SU( 3)f transforms the matter fields in the same 
way as U( 1)pq with an order parameter of 9 = 27r(n/3), which is indeed the Z discrete 
subgroup of C/(1) P q transformation. This means that the Z 3 ® group is identified with the 
center of SU(3){, and thus it is embedded to the continuous S77(3)f group. This implies 
that the degenerate Nf vacua are continuously connected with each other by the SU(3)i 
transformation. 

Here, we comment on the general case in which the number of family is Nf. Let us 
consider a variant of the DFSZ axion model where the interaction term between the Higgs 
fields and SSB breaking field x is given by H u H d x n with a certain integer n. In this case, the 
non-perturbative QCD effects associated with instantons break D(1)pq to Z ^® ■ Here, we 
need to embed Z^ symmetry to the SU(Nf) symmetry to solve the domain wall problem 
by the Lazarides-Shafi mechanism. However, the center of SU(Nf) is Z Nfl so that n has to 
be equal to unity. Then, we should check wether the interaction term H u H d x is consistent 
with the SU(Nf) symmetry. In the case of Nf = 3 it is actually consistent, while in the other 
cases it may not be. This implies that only the case with Nf — 3 can realize the Lazarides- 
Shafi mechanism to solve the domain wall problem. Although the case with Nf = 1 is also no 
domain wall problem, it is disfavoured in light of leptogenesis since there is no CP-violating 
phase. Therefore, our scenario may explain that the number of family is equal to 3. 

3 Constraints and predictions 

In this section, we consider constraints and predictions in our model. Cold axions are gen¬ 
erated around the time of the QCD phase transition from three mechanisms: misalignment 
mechanism [13, 14, 15], emission from cosmic strings [34], and decay of domain walls [35]. 
The resulting axion abundance can be consistent with the observed DM abundance as shown 


( 4 ) 

( 5 ) 

( 6 ) 
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in the next subsection. In Sec. 3.2, we estimate the axion and familon decoupling tempera¬ 
ture and show that they contribute to dark radiation in the Universe [29, 30]. The resulting 
amount of dark radiation can be measured by future CMB observations. Then, in Sec. 3.3, 
we calculate an exotic kaon decay rate that is severely constrained by collider experiments. 
The result is marginally consistent with the present constraint and would be measured by 
future collider experiments. Finally, we check that the seesaw mechanism and leptogenesis 
can be naturally realized in our model. 


3.1 Axion dark matter 


In the case that the U( 1)pq symmetry is spontaneously broken after inflation as we con¬ 
sidered in this paper, cold axions are generated from three mechanisms: the misalignment 
mechanism [13, 14, 15], emission from cosmic strings [34], and decay of domain walls [35]. 3 
The resulting amounts of axions from these mechanisms are the same order with each other. 
The detailed calculations are performed in Ref. [36], though we need some miner corrections 
to apply the results to our model. 

First, let us explain the axion production from the misalignment mechanism [13, 14, 15]. 
At a temperature above the QCD phase transition, the axion field stays at phases that are 
randomly distributed over the horizon scale. As the temperature decreases, it acquires the 
effective mass through the nonperturbative effect [31, 32]. A study based on the interacting 
instanton liquid model implies that the axion mass in a finite temperature can be fitted by 
the following power law formula [37]: 


m a (T) 2 = c T 


A 4 

^QCD 

F„ 2 V A 


T 


'■QCD 


(7) 


where ct = 1-68 x 10 ', a = 6.68, and Aqcd = 400 MeV. This power law formula should 
be truncated when the axion mass reaches the zero temperature value of 


m a ( 0) ~ 


m u m d m-xf-K 
(m u + m d ) 2 F a 


( 8 ) 


When the Hubble parameter becomes to be comparable to the axion mass, the axion field 
starts to oscillate around the minimum of the potential. It occurs at the time around 
t = t osc that is defined by m a (t osc ) = H(t osc ). Since the axion field oscillation is adiabatic, its 
number density is conserved after the oscillation begins. Thus we obtain the axion abundance 

3 A typical momentum of axions generated from topological defects is of the order of the Hubble parameter, 

so that these axions soon become nonrelativistic due to the redshift effect and can be cold DM [36]. 
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generated from the misalignment mechanism as 

nf * (0.07 - 0.09) ( w ^) 119 (j^) . (9) 

where we have substituted 0(1) numerical parameters [36] and h is the Hubble parameter 
in units of 100 km/s/Mpc. Note that the result is independent of initial conditions because 
the initial phase of the axion held is randomly distributed and we should average it over the 
horizon scale to obtain the net axion abundance. 

Next, we consider the axion emission from cosmic strings [34], Since we consider the 
case that the SU(3){ x U(1)pq symmetry is spontaneously broken after inflation, cosmic 
strings form at the time of SSB. The energy density per unit length of cosmic strings /i s tring 
is roughly given by 


/‘'■it.™,, — log (<S s '(), 


( 10 ) 


where S s is the core width of cosmic strings. Cosmic strings emit NG bosons (axions and 
familons) to follow the scaling dependence: 


Pcs - 


/^-string 

I* - ’ 


(ii) 


where p cs is the energy density of cosmic strings. Since cosmic strings lose their energy by 
emitting NG bosons, the energy density of NG bosons is roughly equal to that of cosmic 
strings. Here, the NG bosons consist of familons as well as axions, so that the energy density 
of axions have an additional 0(1) uncertainty in our model. Integrating it from the time of 
PQ phase transition to t osc , we obtain the axion abundance from cosmic strings as 

= ( 0 . 005 - 0 . 2 ) (—^ -^ ( Aqcd 1 , ( 12 ) 

a v ; Vio n Gevy V 4 ooMevy’ V ; 


where we include uncertainties coming from numerical simulations. In addition to that, 
we conservatively include an uncertainty of (1/9 — 1) coming from the efficiency of axion 
emission from cosmic strings, where the factor of 1/9 corresponds to the case that cosmic 
strings emit all types of NG bosons equally. 

Axions are also generated from the decay of domain walls [35]. Although the domain 
wall problem is solved by the Lazarides-Shafi mechanism in our model, short-lived domain 
walls with IVdw = 1 form at the QCD phase transition (t = t osc ). Since the domain walls 
are bounded by cosmic strings, they disappear when the tension of domain wall cr wa n exceeds 
that of cosmic string: 

<Tw.,l(«2) = fe ‘7 (t2) . (13) 

1 2 
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Here, the tension of domain wall is given by cr wall (f) ~ 9.23 m a (t)F%. The energy density of 
domain walls is released as axions when they disappear. Thus we obtain the axion abundance 
from the decay of domain walls as 


Q? w h 2 


(0.004 


0.13) 


10 11 GeV 


1.19 


( Aqcd \ 

\400 MeV J ' 


(14) 


To sum up the above three contributions, we obtain the total amount of cold axions as 

1.19 


ft DM ~ (0.08 - 0.43) 


10 11 GeV 


A 


QCD 


400 MeV 


(15) 


This is consistent with the observed DM abundance ~ 0.12 when the axion decay 

constant is given as 


F a ~ (0.3 - 1.5) x 10 11 GeV, (16) 

where we use Aqcd = 400 MeV. 

3.2 Dark radiation 

While cold axions are generated at the time around the QCD phase transition as we calcu¬ 
lated in the previous subsection, relativistic axions are also produced from thermal plasma 
at high temperature [38]. The relevant and model-independent reactions are scatterings be¬ 
tween gluons and/or quarks in the thermal plasma, whose cross sections are roughly given 
by cr ~ /(87 t 2 F%). A detailed calculation in Ref. [39] reveals that those processes keep 
axions in thermal equilibrium with the thermal plasma until the temperature decreases to 
the following value: 4 * * * * 


T, 


(axion) 


D 


~ 2 x 10 9 GeV 


10 11 GeV 


(17) 


Next, we consider the thermalization of familons. Below the mass scale of y the familons 
interact with the axion through the effective interactions 

C-.H ^ + ■ • •. ( 18 ) 

a f 

4 Reference [40] claims that the hadronic processes such as 7T7t o ttcl are thermal equilibrium below the 

QCD phase transition for F a < 10 12 GeV. However, calculations in Ref. [41] imply that those processes 

are inefficient for F a > 10 s GeV. In this paper, we follow the latter result which allows us to neglect those 

reactions for F a ~ 10 11 GeV, because the estimation by a naive cross section of ov ~ F ~ 2 is inconsistent 

with the former result by more than four order of magnitude. 
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where Ff (= 0(vi)) collectively denotes familon decay constants. Thus the familons are 
produced by the scattering of axions as a + a —> f + f. The scattering cross sections are 
roughly given by cr ~ T 6 / (F*Fj) and those scattering processes are decoupled at temperature 

rgamilon) ~ IQ 9 " 10 GeV. (19) 


If heavy Higgs doublet bosons have a mass scale Mu less than the reheating temperature, 
the Compton-like scattering processes may keep axions and familons in thermal equilibrium 
with the thermal plasma. This process is decoupled at the temperature around T ~ Mu- 
We consider that it is most likely the case, since otherwise the required fine tuning of the 
lightest Higgs boson mass becomes severer. However, whether or not these processes are 
taken into account, relativistic axions and familons are decoupled at a temperature above 
the electroweak phase transition and below the PQ phase transition. Thus, our conclusions 
are not altered whether or not we take into account those contributions. 

As we show in Sec. 3.4 the thermal leptogenesis requires T rh > 3 x 10 9 GeV, so that 
axions and familons are generated from the thermal plasma after reheating. Then they are 
decoupled well before the electroweak phase transition and contribute to dark radiation in 
the subsequent cosmological history [29, 30]. The amount of dark radiation is conventionally 
expressed by the effective neutrino number N eS defined by 


Prel 


P'1 


1 + N t 


eff: 



( 20 ) 


The deviation from the standard value can be calculated from 


AV eff - ^ ~ 


4 f 9*{T' d ) 


7 V 43/4 


-4/3 


( 21 ) 


where the summation is taken for all relativistic particles i decoupled at T l D . Using g*(Tu) = 
106.75, we obtain A N^ lon ~ 0.027 for the axion contribution to the effective neutrino 
number. Taking into account familons as well as axion, we predict the amount of dark 
radiation as 


A N eS ~ 0.24. (22) 

An extra background of relativistic particles affects the CMB anisotropies by the early 
integrated Sachs-Wolfe effect. The abundances of light elements are also affected by dark 
radiation through the effect on the expansion rate at the BBN epoch. The Planck data 
combined with the observation of helium abundance estimated by Ref. [42] is consistent 





with the standard value N e g = 3.046 but still allows additional radiation component with 
the constraints [43] 


N eS = 2.99 ± 0.39 (95 %CL). 


(23) 


The measurement of deuterium abundance [44] is consistent with this result. On the other 
hand, Ref. [45] have reported a larger central value from the measurement of helium abun¬ 
dance: 


N eS = 3.58 ± 0.40 (95.4%CL). 


(24) 


However, it is under discussion whether the uncertainty accurately reflects systematic errors. 
Our result of A N e g ~ 0.24 is consistent with those observations. 

Future observations of the CMB fluctuations will measure the effective neutrino number 
with precisions of 


A N eS = 0.0156 (lcr), 


(25) 


by the ground-based Stage-IV CMB polarization experiment CMB-S4 [46] (see also Ref. [47]). 
We conclude that familons and axions will be indirectly observed as dark radiation in the 
near future. The CMB-S4 experiment would indirectly observe axions even if dark radiation 
consists of only axions (AN e g = 0.027), which is the case of the KSVZ axion model [8], for 
example. Let us emphasize that the measurement of the amount of dark radiation will give 
us a probe of the number of NG bosons [29, 30]. 

3.3 Collider signatures 

In this subsection, we consider collider signatures predicted in our model. In general, NG 
bosons f a associated with symmetry breaking interacts with matter fields via 



(26) 


where j afl is the broken symmetry current and Fj is the decay constant. In light of the 
collider experiments, the most relevant interaction is given as 5 


~^d fl f 2 sYd + h.c ., 
b h 


(27) 


5 Although the familon interactions also lead to flavour changing lcptonic decays such as n — > ef and 
/i —> ey/, these decay modes pose weaker lower bound on Ff than the K + decay mode does [48]. We should 
also note that it is very challenging to improve the limits on these leptonic decay processes with the current 
experimental facilities [49]. 
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where Fj 2 is given by \f2 |t>i — u 2 |. 6 The flavour changing interaction of Eq. (27) leads to an 
exotic kaon decay process K + — > 7 r + /. The reaction rate is calculated as [48] 

r(A' + ^ I r + /) = 4-^4/J 3 |C 1 (0)| 2 , (28) 

h 

where (3 = 1 — ml/m 2 K . The form factor (ti + (p') |s 7 M d| K + (p)) = Fi(q 2 )(p + p'Y has a 
normalization Aj(0) = 1 in the limit of exact flavour SU(3) symmetry. Collider experiments 
pose upper bound on the exotic kaon decay rate as [50]' 

Br(A" + -A 7 r + f) < 1.3 x 10” 10 (90%CX). (29) 

This leads to a lower bound on the decay constant: 

F h >5.3 x 10 11 GeV. (30) 

This is marginally consistent with the decay constant F a determined by the cold axion 
abundance (see Eq. (16)) because F h (=V2\ Vl — 772 !) can be larger than F a given by Eq. (3) 
when one of v l is smaller than the others. 

The recent start of the NA62 experiment at CERN is expected to measure the branching 
ratio of the process K + —> n + vv within an accuracy of 10% compared to the SM prediction 
Br(A' + —» 7r + z/i7) SM = (9.11 ± 0.72) x 10~ n [52] by 2018 [53]. This implies that the NA62 
experiment will observe a deviation from the SM prediction when the decay constant is just 
above the present lower bound. 


3.4 Neutrino mass and leptogenesis 


In this subsection, we introduce right-handed neutrinos v c to explain the small masses of left- 
handed neutrinos by the seesaw mechanism [10, 54], Since they are charged under SU( 3)f 
as shown in Table 1, their interaction terms are written as 


C 




(31) 


where y and y' are Yukawa coupling constants. After the SSB of SU( 3)f x U( 1)pq, the 
Lagrangian is reduced to be 


C 


--M£M-y ii L i H u rf + h.c., 


(32) 


6 Here we implicitly assume that the left-handed flavour eigenstates (c?l,sl) interact with familons in 

the same way as the right-handed ones (d_R, sr). In general, this is not the case and we should include the 

flavour mixing effect. In this paper, we neglect the flavour mixing effect for simplicity. 

'Although Ref. [51] have reported the slightly severer upper bound as Br(/\ + —> tt + /) < 4.5 x 10 —11 , we 

conservatively quote the upper bound derived in the recent paper [50]. If we take the former upper bound, 

the lower bound on Ff 2 is calculated as Ff 2 ^ 9.0 x 10 11 GeV. 
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where = y (x y ). Here, we have integrated out heavy Higgs bosons and have written 
the light Higgs H l J as y'H l J = y l] H u . Integrating out u c , we obtain the mass matrix of 
left-handed neutrinos via the seesaw mechanism: 


m u = —m D M u ^m^), (33) 

where (m^)*- 7 = y l i ( H u ). The observed neutrino mass deferences ~ (7.53 ± 0.18) x 

10~ 5 eV 2 and A m\ 2 — (2.44 ± 0.06) x 10 -3 eV 2 (normal mass hierarchy) [55] imply that 
y* = C>(10- 3 ) for M% = C>(10 10 ) GeV. 

The baryon asymmetry of the Universe can also be explained by the thermal lepto- 
genesis [11], We assume that the reheating temperature is larger than the mass of the 
lightest right-handed neutrino M u c. Assuming hierarchical right-handed neutrino masses, 
M v c M u c , we can calculate the amount of baryon asymmetry via 


v n b 
Yh = — ~ 


24 + 4n 


H 


enY eq , 


(34) 


*• 66 + 13 uh 

where e is the CP-asymmetry parameter in decays, and the fractional factor, depending on 
the number of Higgs doublet n# (= 2 in the DFSZ model), comes from the sphaleron effect. 
The equilibrium abundance is given by Y eq = 135([(3)/(47r 4 g*), where g M is the effective 
number of spin-degrees of freedom in the thermal plasma. The efficiency factor k has to be 
included to take into account the washout effect and the correction of the number density 
of vf with respect to the equilibrium value U eq . 

The CP-asymmetry parameter is calculated as [56] 


e ~ 


< 


3 -M„ 


m 




167T 


(Hu) 


2 ) 


(35) 

(36) 


where the last inequality is known as the Hamaguchi-Murayama-Yanagida bound [57] (see 
also Ref. [58]). It is convenient to define the effective neutrino mass fh\ by 


mi = 




(37) 


The lightest left-handed neutrino mass m v c satisfies the inequality of < rhi. The effi¬ 
ciency factor k has been obtained from a detailed calculation in Ref. [59] and has been fitted 
by 


K 


-1 


3.3 x 11) ° ev 

m i 


/ mi \ 

V0.55 x 10- 3 eV ) 


(38) 
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The right-handed neutrino cannot be efficiently produced from the thermal plasma for 
the case of fh\ <C 1CV 3 eV, while washout effects are efficient to dilute lepton asymmetry for 
the case of rhi 10~ 3 eV. These are the reason that the efficiency factor has a maximal 
value 0.2 around rh\ ~ 10” 3 eV. 

The observations of CMB fluctuations imply that the baryon-to-photon ratio is given 
as [55] 

— = (6.05 ± 0.07) x 1(T 10 , (39) 

n 7 

where the number density of photon n 7 is related to the entropy density as s = 7.04n 7 . Using 
Eqs. (34), (36), and (38), we can explain the amount of baryon asymmetry by the thermal 
leptogenesis for the case of 7 r h ~ M u c ^3x 10 9 GeV (see, for a review, Ref. [23]). Note that 
the masses of the right-handed neutrinos come from the VEV of y, which is related to the 
axion decay constant F a . The result of M v c > 3 x 10 9 GeV is consistent with the result of 
F a = (P(10 n ) GeV obtained in Sec. 3.1. The reheating temperature of Trh ~ 3 x 10 9 GeV 
is also consistent with our assumption of the thermalized axion and familons, since they are 
decoupled at temperature around T 10 9 GeV. 

4 Conclusion 

We have proposed a QCD axion model with a SU( 3)f flavour symmetry to solve the domain 
wall problem by the Lazarides-Shafi mechanism. The unwanted discrete PQ symmetry is 
embedded to the continuous SU( 3)f symmetry to connect the vacua. Since the SU(3){ 
symmetry is anomalous, it should be global symmetry and predicts eight NG bosons called 
familons in the low energy effective theory. We find that familons as well as axions are 
thermalized and then decoupled after the SSB of SU (3)f x[7(1 )pq symmetry. They contribute 
to dark radiation in the subsequent cosmological history and the resulting amount of dark 
radiation will be detected by future observations of CMB fluctuations. Our model also 
predicts a sizable exotic kaon decay rate, which is marginally consistent with the current 
upper bound. We expect that collider experiments would observe excess of kaon decay signals 
such as K + —>■ 7 r + f in the near future. 
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